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Abstract 

We further develop the approach recently used to construct an analytic ghost-free 
model for the QCD running coupling based on the requirement of the Q 2 -analyticity and 
apply it to the process of e + e~ annihilation into hadrons to study the renormalization 
scheme dependence of the R(s) cross-section ratio. 

By transforming the relevant QCD corrections up to the three-loop level into the "an- 
alytized" form we show that the Raa(s) expression thus obtained is remarkably stable (as 
compared to the conventional perturbative approach) with respect to the renormalization 
scheme dependence for the whole low-energy region. 



1. The analytic model for the QCD running coupling has been devised in Ref. [l] on 
the basis of the Kallen-Lehmann analyticity (that is causality). This expression a an ((5), by 
construction, has no unphysical singularities, like a ghost pole, in the complex Q 2 -plane and 
its perturbation expansion precisely coincides with the usual perturbation one. The simplest 
version (discussed in detail in Ref. @) has no extra parameters. Its infrared limit as Q 2 — > is 
finite and independent of the QCD scale parameter A. The numerical value a an (0) ~ 1.4 turns 
out to be remarkably stable with respect to higher loop corrections. 

The analytic running coupling, by prescription, is defined via the spectral integral 
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with the density p(a) calculated "perturbatively" as an imaginary part of the usual (RG- 
summed) running coupling. This spectral representation provides us with correct analytic 
properties of the running coupling a &n (Q 2 ): it is an analytic function in the complex Q 2 -plane 
with a cut along the negative part of the real axis. 

The one-loop expression for the analytic running coupling thus defined can be presented 
explicitly 
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with (3\ = (11 — 2//3)/4 and /, the number of active quarks. The ghost pole "is removed" by 
the second term, the term that is nonanalytic in the coupling constant a M = (p\ In /x 2 /A 2 ) -1 , 
i.e., has a nonperturbative nature. 

Two- and three-loop expressions for a an (Q) discussed in our previous publications have a 
more complicated structure and can explicitly be presented only after some approximation^ 
Their analysis has revealed an important feature - reasonable stability of the analytic coupling 
behaviour in IR with respect to higher loop corrections. A very small difference between two- 
and three-loop curves (see Figure at Ref. @) gives, in general, hope of stability with respect to 
the scheme dependence. 

In this letter, we concentrate on the problem of the renormalization scheme (RS) depen- 
dence. Along with the discussion in Ref. |3|], we are going to apply the idea of analytic approach 
(AA) not only to QFT objects, like a running coupling, but directly to observables. Here, we 
consider R(s), the well known cross-section ratio for the process e + e~ — > hadrons. This ex- 
ample is physically interesting for the RS stability issue (see, e.g., discussion in Ref. Q and 
references therein). At the same time, it is instructive from the theoretical point of view, as 
the time-like region is involved. 

2. The QCD correction r(s) to the physical (and, therefore, RS-invariant) quantity^ R(s) 
is usually written down 

R(s) = 2[l + r(s)]; r(s) = a(s) [l + r x a(s) + r 2 h 2 (s)] (3) 

in terms of a(s), "the running coupling in the time-like region" which, by assumption, is defined 
as a mirror image a(— Q 2 ) = a(Q 2 ). This expansion is parallel to the analogous one for the 

^or the two-loop case one can use expression (|^) with the substitution 

Q 2 /A 2 -> exp 



In Q 2 /A 2 + /3 2 /(/?i) 2 ^ln 2 Q 2 /A 2 + 4^ 2 



see Ref. 



2 Here and below all numerical coefficients are given for three active quarks that is adequate to the "low-Q 
physics" we are interested in. 
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RS-invariant Adler function 

D{Q 2 ) = Q 2 ^^ = 2[l + d(Q 2 )], 



dQ 2 



d(Q 2 ) = a(Q 2 ) [1 + d ia (Q 2 ) + d 2 a 2 (Q 2 )\ , (4) 
where Tl{Q 2 ) is the vector current correlator connected with R(s) by the relation 

Here, r\ = d\ but the second (and higher) coefficient differs r 2 — d 2 = — (/?i7r) 2 /3 = — 16.65 
by the so-called "7r 2 -term" which is the "trace" of an analytic continuation from the space-like 
to the time-like region. Taking into account that d(Q 2 ) is an analytic function in the complex 
Q 2 -plane with a cut along the negative part of the real axis, the analytic continuation procedure 
can be described by the "reverse relation" 

/ \ 1 f s+i£ da 

r{s)=-— —d(-a) } (6) 
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where the contour is in the region of analyticity of d(Q 2 ). 

However, this straightforward combination of analytic continuation with the "mirror" RG 
summation preserves unphysical singularities, like a ghost pole in the one-loop term and un- 
physical cuts for higher loops. Instead of it we shall apply the AA. 

To compare the results obtained in various RS's, one uses the 'cancellation index criterion' 
proposed by P. Raczka |J. According to it, a set of "natural" RS's can be introduced which 
includes schemes for which the degree of cancellation between different terms in the second 
RS-invariant || 

p2 = r 2 + c 2 - r\ - ci t x = p{ci, r f ) (7) 
is "not too large" . The degree of cancellation is measured by the 'cancellation index' 

C R = T—r ( |r 2 | + |c 2 | + r\ + ci M) = C(c h n) , (8) 
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and the usual quantitative criterion is || 

C R <2. (9) 

Note here, in order to compare various RS's, one can also use, instead of index Cr, Eq. (g), 
the "Euclidean index" Ce = C(ci, dj) expressed via the Adler function coefficients d& that are 
more primary than as they emerge directly from the Feynman diagram calculation. This 
change Cr — ► Ce with appropriate redefinition of the second RS-invariant p 2 — > p(cj, dj) (see 
Ref. 0) will influence the classification of "well-behaved schemes". We shall not discuss here 
this point in detail. 

3. To implement the AA, we will base upon Eq. (j5|) for QCD correction to the Adler 
function and exploit the analytic property of d(Q 2 ) expressed in the form 

d(Q 2 ) = - r-^-^(a). (10) 
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It is essential that the Adler function is defined in the Euclidean region where the renormaliza- 
tion group (RG) method can be applied directly. 

On the other hand, using the reverse relation for r(s), Eq. @, we have (see Ref. f|) 

r AA (s) = - r d —?\a). (11) 



TC Js a 

Here, we introduce an effective spectral density g eS (cr) which corresponds to the disconti- 
nuity of d(Q 2 ) and can be presented as 

g cS {a) = p {a) + d 1 6l (a) + d 2 g 2 {a) , (12) 

where the first term is just the density for the running coupling and Qk{cr) are related to higher 
powers ~ a k+1 (Q 2 ) in the expansion with the coefficients || d?f s = 1-64 and g^ s = 6.371 (the 
latter corresponds to r^ s = —10.284). 

The coefficients and dk in Eqs. (fl2|) and (f|) as well as c 2 in the running coupling a{Q 2 ) 
depend on the RS choice. In the ftfunction 

Q 2 J^ = -A a 2 (l + c x a + c 2 a 2 ) ; a = a(Q 2 ) , (13) 

we have ft = 9/4 = 2.25, c x = 1.778 and, e.g., cf s = 4.471. From Eq. ([TJ), taking into account 
the first RS-invariant p\ = ft \ia(Q 2 / A 2 ) — d\ and using a scale parameter A^jg in the MS scheme 
as a reference one, one finds the relation for the running coupling 

ft In = <^ - d\ + $(a*(Q 2 ), 4) (14) 
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$(a c ) — * + c In + c r ^ X (15) 

~ a (1 + ci a) jo (1 + cix) (1 + c\x + C2X 2 ^ 



and i, the RS index. 

To calculate the QCD correction r(s) within the conventional approach, one takes Eq. ( II ) 
with Q 2 — > s to define the "running coupling" a(s) in the time-like domain and substitute it 
into the expression for r(s), Eq. (||). As it has been noted above, this procedure encounters 
difficulties with analytic properties of the usual running coupling. Within the AA, we arrive 
at the same quantity, taa(s), from Eq. (0). Here, the effective spectral function g eS (a) is 
defined as an imaginary part of the expression (|J), with the function a(Q 2 ) being replaced by 
the complex function a(— cr) which has to be found numerically as a solution of Eq. ( [[4]) with 
the replacement \n(Q 2 /A 2 j^) — > ln(cr/A|^) — m. 

To discuss the RS-dependence, we consider two examples of "well-behaved" RS's: the 
scheme A with the parameters r x = —3.2 and = (the so-called 't Hooft scheme) and 

MS scheme with rj MS ^ = 1.64, c 2 MS ^ = 4.47. Both the schemes have a sufficiently small value 
of the cancellation index Cr ~ 2, which is close to that for the so-called optimal RS based on 
the principle of minimal sensitivity (see the discussion of this subject in Ref. [|i~T| ). 

3 The explicit expression for <I>(a, c?) can be found, for example, in Ref. JlCf] . 
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Figure 1: Plot of the QCD correction r(s) calculated in the cases of perturbation theory (PT) 
and the analytic approach (A A) in two different RS's with approximately the same cancellation 
index C R ~ 2: A and MS. 

In Fig. [I], we plot the QCD correction r(s) as a function of v^M-ms for these two schemes 



in usual treatment, as it was considered, e.g., in Refs. fj, [l^, [H| and within the AA. One can see 
that the analytically improved result for R(s) obeys a stable behaviour for the whole interval 
of energies being practically scheme independent. 

Note also that the values of the parameter A in the AA and in perturbation theory are not 
equal to each other due to nonperturbative contribution to the AA running coupling for the 
same value of the coupling constant at some reference point. The value of Aaa is larger than 
Apt by a factor of about two 0. Hence, our AA curves and the MS one, being considered as 
functions of ^/s, will be rather close to each other at s > 2 GeV. 

4. In our previous paper 0, we have demonstrated that the analytically improved running 
coupling being a smooth function in the IR turns out to be remarkably stable there with respect 
to higher loop corrections. 

In this note, we have applied the same idea of the analytical approach to study the problem 
of theoretical ambiguity coming from the RS dependence for a physical quantity, R(s), the cross- 
section ratio of the process e + e~ to hadrons. In the next-to-next-to-leading order the Adler 
function with the correct analytic properties has been constructed. This fact allowed us to find 
the corresponding "analytized" expression for the QCD correction r(s) to R(s). In contrast 
to the conventional perturbative expansion this procedure does not encounter any problems 
and can be performed in a self-consistent way. We have found that our AA reduces the RS 
dependence drastically: Raa( s ) thus obtained turns out to be practically RS independent for 
the whole energy interval. 
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